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The relation between three-dimensional lattice structure and magnetism in correlated electron 
systems is explored for face centered cubic (FCC), body centered cubic (BCC), and simple cubic 
(SC) lattices. In particular, we question which lattice structure has the strongest tendency 
toward the ferromagnetism or antiferromagnetism. We employ the Hubbard model to calculate 
the spin susceptibility and the single-particle spectrum with the fluctuation-exchange (FLEX) 
and the two-particle self-consistent (TPSC) approximations in the weak coupling regime. We 
have shown that (i) ferromagnetic spin fluctuations become dominant when the Fermi level 
lies around a sharp peak (divergence) in the density of states (D(E)) near the bottom of the 
band, which occurs for FCC with/without next nearest neighbor hoppings (t 1 ) or BCC with an 
appropriate value of t 1 . Among the cases studied, the ferromagnetic fluctuation is found to be 
the strongest for FCC with a finite t'. (ii) When the peak in D(E) resides around the band 
center as in bipartite SC or BCC, antiferromagnetic fluctuations become dominant when the 
band is close to the half- filling, with the fluctuation being much stronger in BCC. 

KEYWORDS: three-dimensional Hubbard model, ferromagnetism, antiferromagnetism, fluctuation exchange ap- 
proximation, two-particle self-consistent approximation 



§1. introduction 

The metallic ferromagnetism in repulsively interact- 
ing electron systems is among central problems in solid 
state physics, but is still some way from a complete un- 
derstanding despite a long history of investigations. Af- 
ter the seminal papers by Kanamori,El by Hubbardp 
and by Gutzwillerff appeared back in the 1960's, the 
single-band Hubbard model has been studied intensively 
to understand metallic ferromagnetism. Although it is a 
simplest one conceivable, the model can still harbor in- 
sights into the fundamental mechanism for an itinerant 
ferromagnetism in correlated electron systems. 

Historically, a guiding principle for realizing ferromag- 
netic states was provided by Stoner's mean field theory, 
which suggests that the ferromagnetic state should be 
favored for large enough interaction and/or density of 
states (DOS) at the Fermi level. KanamoriEP elaborated 
this argument with the T-matrix approximation to show 
that the ferromagnetic states are indeed favored when 
a peak in the DOS around the Fermi level is embed- 
ded around the bottom of the wide band. However, for 
general band filling, the extent to which the T-matrix 
approximation is valid is not clear. 

As a different avenue along this line, rigorous results 
were recently obtained for ferxomagnetism on specific lat- 
tices having flat bands. Liebtf has shown that, when we 
have a bipartite lattice with different numbers of sub- 
lattice sites, the ground state at half filling should be 
ferrimagnetic for arbitrary magnitudes of the repulsion, 
U. Different numbers of sublattice sites imply that there 
is a flat band(s) in the one-electron band structure, or 
a delta-function density of states. Mielke and TasakitJ 



then proved the existence of ferromagnetic ground states 
in systems with flat band for more general classes of lat- 
tices. 

Going back to ordinary lattices, one can then look for 
divergences in the DOS to realize a ferromagnetic state 
in the Hubbard model. For three-dimensional lattices, 
the divergence occurs for the face centered cubic (FCC) 
lattice (where the DOS diverges at the band bottom, 
or above the bottom if we include next-nearest neigh- 
bor hopping, t') or for the body centered cubic (BCC) 
lattice (where the DOS diverges at the center due to 
the electron-hole symmetry, or away from the center 
if we include an appropriate value of t'). Our ques- 
tion here is 'among the lattice structures having sharply 
peaked (divergent) DOS's which one has the strongest 
tendency toward ferromagnetism?'. While the fully po- 
larized ground state is expected to be favored when the 
DOS diverges near the bottom of the band, the tendency 
toward the ferromagnetic instabilities at finite tempera- 
tures is a non-trivial problem, because Tc of a system 
whose ground state is fully polarized is not necessarily 
higher than that of a system whose ground state is par- 
tially polarized. 

As the background for the problem we can first reca- 
pitulate how the link between the lattice structure and 
the ferromagnetism has been understood in terms of the 
Hubbard models' The link may be tested accurately 
in one-dimcnsional(lD) systems, since the density ma- 
trix renoraialization group (DMRG) method developed 
by WhitaJ' is applicable in ID. Recently, a saturated 
ferromagnetic ground state has been found in the Hub- 
bard jmodfil on various lattices. These include a trestle 
chaiBP'Q&3> which may be regarded as a generalized flat 
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band model in 1A0* ladder systemstBli^ and nearly 
flat-band models xS 

In the opposite limit of infinite dimensions, exact 
calculation by means of dynamical mean field theory 
(DMFT) becomes feasible, since the self energy becomes 
site-diagonal or k-independent. Using finite-tempexature 
quantum Monte Carlo (QMC) techniques, UlmkeEiJ has 
solved the DMFT equations to discuss the stability of 
metallic ferromagnetism for .ECC-type lattices in infinite 
dimensions. Wegner et al$3 obtained similar results 
with an iterated perturbation theory instead of QMC. 
If one compares these results with, that for the hyper- 
cubic lattice by Obermeier et al.,U$ FCC is seen to be 
favorable for ferromagnetism. 

For two-dimensional (2D) or three-dimensional (3D) 
systems, which are of realistic interests, the situation is 
less conclusive, since we cannot apply thiiLL= oo result 
straightforwardly. For 2D Hlubina et alS3x3 suggested 
a possibility of the ferromagnetic ground state when the 
Fermi level is located near the van Hove singularity in 
the t-t' Hubbard model for low densities. For 3D there 
are many studies with a variety of approximate methods. 
Hanisch et al& studied the instability of Nagaoka's sat- 
urated ferromagnetic ground stateEd' for various lattices 
with a variational method. The DMFT, which has origi- 
nally been developed for d — oo, has also been applied to 
3D by Ulmkdj by plugging the DOS of the FCC lattice 
into the DMFT equations to obtain the phase diagram 
for FCC. However, the validity of the approximation, i.e., 
whether we can neglect the k-dependence of thp_self en- 
ergy in 3D, is quite an open question. NoltingEa) et al. 
also studied the stability of the ferromagnetic state for 
the FCC Hubbard model with the spectral density ap- 
proach, which JHay give reasonable results in the strong 
coupling limit ,£5 but the obtained Tc turns out to be 
several times larger than those estimated in the DMFT. 

Since the above approaches involve one approximation 
or another (such as the strong-coupling limit, infinite di- 
mension, etc), further investigations with different ap- 
proaches should shed a new light on the problem. The 
purpose of the present paper is to investigate the mag- 
netic properties of the 3D Hubbard model from the weak- 
coupling side, where the fluctuation-exchange (FLEX) 
and the two-particle self-consistent (TPSC) approxima- 
tions are used to explore ferromagnetism for the first 
time. With these methods we look into FCC, BCC, and 
simple cubic (SC) lattices with or without next-nearest 
neighbor hoppings. r~ir~i 

The FLEX, introduced by Bickers et a/.,E§EZP takes 
a set of skeleton diagrams for Luttinger-Ward func- 
tional to generate a (fc-dependent) self energy that is 
computed self^cnnsistently based on the idea of Baym 
and Kadanoff.C-jEj) Hence the FLEX is a kind of self- 
consistent perturbation approximation with respect to 
U. The method has^mL-applied to the analysis of the 
2D Hubbard model BEZ1e#§) in the context of high-T c 



cuprates, superconducting ladder systemsEf or organic 
superconductors.EirEj-' The FLEX breaks down in the 
intermediate to strong coupling regimes where the ver- 
tex corrections become significant, as seen from the fact 
that the pseudo-gap or the Hubbard bands are not repro- 



duced, c3> but the method should give reasonable results 
at least in the weak-coupling regime, namely as far as 
the fluctuations do not grow too much. 

In Kanamori's theory jP vertex correction is taken into 
account with the T-matrix approximation, which is only 
valid in the dilute limit. Here we have opted for the 
FLEX approximation, since the theory, with no restric- 
tions on the electron density, can discuss the density 
dependence of the spin fluctuation, although the ver- 
tex corrections are neglected with only the self-energy 
corrections are taken iiito account there. In this context, 
Moriya and TakahashEJ developed a theory which takes 
account of both the vertex corrections and the selfpenergy 
corrections. Along this line, Usami and MoriyaH dis- 
cussed the ferromagnetic metals such as Fe, Co, Ni with 
a functional integral method. Similar xtr more elabo- 
rate calculations were done by HasegawatS' and by Hub- 
bard S3 However, these are done at the cost of introduc- 
ing phenomenological parameters, while the FLEX can 
be implemented from a microscopic level (e.g., the Hub- 
bard model) without introducing any phenomenological 
parameters. 

The TPSC method, on the other hand, takes vertex 
correction into account in a self-consistent manner at the 
two-body level by assuming a parameterized form for the 
Luttinger-Ward functional, where renormalized interac- 
tions U sp and £/ch are introduced. In term of these 'effec- 
tive [/' we assume the same functional forms (see eqns 
(19,20) below) as in the RPA for charge/spin susceptibil- 
ities. These parameters are determined so as to satisfy 
the sum rules or the constraints on a two-particle prop- 
erty, e.g., the double occupancy {n^n\). Eor the 2D Hub- 
bard model on square lattice, Vilk et al& succeeded in 
reproducing the pseudo-gap and the upper/lower Hub- 
bard bands, which indicates that the TPSC approach 
may._be valid in the intermediate coupling regime as 
well.E3 However, a numerical difficulty can arise in the 
calculation of the chemical potential as we shall see later 
(§IIIA), so we concentrate on the FCC lattice with next- 
nearest neighbor hoppings, for which the ferromagnetic 
spin fluctuation turns out to be the strongest among the 
cases studied in FLEX, and is of interest. 

Both FLEX and TPSC methods will become invalid in 
the strong-coupling regime, while strong repulsion is of- 
ten thought to be required for the realizatiotLof ferromag- 
netic ground states as in Nagaoka's caseO Neverthe- 
less, we expect that the tendency toward a ferromagnetic 
phase transition in the weak-coupling approach may be 
captured even when the transition itself occurs outside 
the regime. Here we have numerically checked to which 
value of U the method may be relied upon. Another pur- 
pose of the present paper is to discuss in which lattice 
structure the aniiferromagnetism is most favored. Anti- 
ferromagnetism, contrary to ferromagnetism, is tractable 
with approaches in the weak-coupling regime in general. 

The present paper is then organized as follows. In 
section II, after introducing the model, we review the 
formalism of FLEX and TPSC. Section III is devoted to 
the results and discussions, while section IV summarizes 
the results. 
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§2. Formulation 

2.1 Model Hamiltonian 

We consider the conventional single-band Hubbard 
model, 

H = % 4r C J* + U J2 n iT n 4> (!) 

where c\ G creates an electron at the i-th site with spin 
a, m a = c\ a Cia is the number operator. We take — 
t for nearest neighbors and tij = t' for second-nearest 
neighbors. The energy dispersions for FCC, BCC, and 
SC lattices are given as 

3 

e FCC (k) =4t^cos(A*)cos(fc 3 -) + 2t'^cos(2fci), (2) 

i<j i—1 

3 

e BCC (k) = 8tcos(fci)cos(fc 2 )cos(fc 3 ) + 2t' cos(2fc l jj3) 



e sc (k) = 2t^]cos(fc l 



(4) 



respectively, where fci = k x , &2 = k y , k% = k z . To 
facilitate the calculation here we take a cubic Brillouin 
zone (— 7r < ki < n) by considering two(four) equivalent, 
interpenetrating FCC (BCC) lattices. Hereafter, we set 
t = 1. 

2.2 The FLEX approximation 

We first review the formalism of the FLEX approxi- 
mation. We start from the Luttinger-Ward functional $ 
which is defined as £1 = $ + T ^ [In G - EG] , where SI is 
the thermodynamic potential, G is the dressed Green's 
function and E is the self energy. The Feynman diagrams 
for <!> considered in the FLEX approximation are shown 
in Fig. 1. 

Fig. 1. The Luttinger-Ward functional considered in the FLEX 
approximation. Full lines represent the dressed Green's function 
while broken lines U . The summation is over n, the total number 
of interaction lines. 



The self energy can be obtained with a functional 
derivative as 



E(fe) = 6<f>/5G 



N ■ 
i 

+G(k - q)V {ph) (q) + G(-k + q)V ipp) (q)} .(5) 

This consists of the contribution from the bubble dia- 
grams, 

v^(q)^u 2 Xpi M, (6) 

the contribution from the particle-hole channel, 

1 



V^{q) = \u\M 



l + Ux P h{q) 



1 



+^ 2 x P h(<z) 



i -Uxph(q) 



- 1 



(7) 
(8) 



and the contribution from the particle-particle channel 

1 



V^(q) = -U 2 Xpp (q) 



where 



Xph(g) 



l + U XpP (q) 
I^G(fc + g)G(fc), 



1 



(9) 

(10) 

(11) 
(k,kj„), 



Here we have denoted q = (q, ie u ) and k 
e v = 2nvT is the Matsubara frequency for bosons while 
uj n = (2n— 1)ttT for fermions, and N is the total number 
of sites. 

The Dyson equation is written as 



G(k,w„) 



G°(k, w „) 



E(k, w n ), 



where 



G°(k, w „) = 



(12) 



(13) 



is the bare Green's function with e° being the energy for 
the free electrons. 

We solve the equations (5) ~ (|l^ ) by setting the chem- 
ical potential fj, so as to fix the density of electrons. The 
Green's function is computed by iteration until a conver- 
gence condition, 



\G {r) (k) - G (r - 1} (fc)|/|G (r) (fc)| < 1.0 x 10" 



(14) 



for G^(fe), the Green's function at the r-th iteration, 
is attained for all the values of k. We have checked the 
convergence by taking (i) 16 3 sites for which 1024 Mat- 
subara frequencies (-1023ttT < e„ < 1023ttT) or 2048 
frequencies are considered, and (ii) a larger 32 3 sites with 
1024 Matsubara frequencies. 

We define the RPA spin susceptibility as 

There is no vertex correction in this expression. As 
for the static magnetic susceptibility, x(0) f can be cal- 
culated, as pointed out by McQueen et a/.,EP from the 
slope Tr [a z G] versus small homogeneous field h (typically 
0.005t ~ O.Oli), where a z is the z-component of the Pauli 
matrix. This is equivalent to solving an integral equation 
for the response function with irreducible vertices com- 
prising particle-particle and particle-hole bubble chains 
along with the Aslamazov-Larkin type diagrams. 

The integrated spectral function, or the DOS, is given 

by 

p{uj) = ^2 ImG(k, uj + iS). 

k 

We obtain G(k, u) by a numerical analytic continuation 
from the imaginary-axis data, G(k, iw n ), with the Pade 
approximation HP 
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2.3 The TPSC approximation 

Next we review the TPSC method following the ar- 
gument by Vilk and TremblayEJ This methods starts 
from assuming that the Luttinger-Ward functional $ can 
be parameterized by two variables, T^-a- and T aa , as 



$ = 



a 

+G CT (0 



r CTCT G CT (o+) 
rv ff G_ CT (o+) 



with 



G ff (0 H 



AT ^ 

fe 



exp{ikO+)G <J (k) 



(15) 



(16) 



If we introduce spin and charge irreducible vertices as 



U sp 



(17) 
(18) 



the charge/spin susceptibility must satisfy the sum 
rules, 



T 
N 



-E- 

N ^ 1 

i 



2x P h(g) 



^chXphte) 
= n + 2(n T n i ) - n 2 , 



n J2 Xs p^ ~ AT E i 

9 <? 



2x P h(q) 



t / s P Xph(g) 



= n - 2(n T n x ). 



(19) 



(20) 



The RPA-like form of the susceptibilities can be read- 
ily derived if we assume that the irreducible four-point 
vertices T's are local in time and space. Since the self 
energy corresponding to the trial Luttinger-Ward func- 
tional is constant, the irreducible susceptibility x p h co- 
incides with the non-interacting form. As a first step 
in the approximation for the chemical potential /i, we 
choose the same value as that of the non-interacting sys- 
tem, flQ. 

The equations (19) and (20) determine U sp and U c h as 
a function of the double occupancy, (n-|-n^). If we follow 
the idea of the local-field-approximation for the electron 
gas due to Singwi et al.juj' we can assume a relation, 



<m(o) 



(21) 
(22) 



from which U sp or U c h can be determined self- 
consistently. 

In the next level of the approximation, the self energy 
is assumed to be 

X [Us P Xs P (q) + U chXcb (q)} G<?\k + #3) 



for which Green's function is given as 

1 



G«(fc) = 



(24) 



where /ii is obtained by the condition, 

n=^GW(fc)exp(-ic n 0-). 

k 

An important advantage of this expression for the self 
energy is that it satisfies the relation, TrE^^G^ -' = 
2U{n-[ni). We can then use the extent to which this 
relation is numerically fulfilled as an internal accuracy 
check. Using G (1 \ Vilk et alE$ succeeded to reproduce 
the pseudogaps and the Hubbard bands in the 2D Hub- 
bard model on square lattice. 
When the Luttinger sum rule, 

1/A^ #(-e k + Ml - E«) = r V , 

is satisfied, we can take the approximation [i\ — 
S(kp, 0) = Vilk et al.tB suggests that, if one consid- 
ers a region near the half-filling or in the strong-coupling 
regime, one may improve the approximation by using 
/ii — E(kp, 0) in place of ^ in calculating the irreducible 
susceptibility, Xph- 

§3. Numerical results 

3.1 FCC lattice 

Now we turn to the present numerical results. We first 
focus on the FCC lattice. We study the following two 
cases; (i) a finite t 1 = 0.5 with n = 0.2,0.6,1.0 and (ii) 
nearest-neighbor hopping only (t' — 0) for the same set 
of band fillings. The DOS for the non-interacting case 
is displayed in Fig. 2(a), which is highly asymmetric 
and has a peak around the bottom of the band. This 
is exactly the situation studied by KanamoriEP with the 
T-matrix approximation: a dilute case with a large DOS 
at the Fermi level. 



Fig. 2. The density of non-interacting states for (a) FCC lattice 
with t! = (solid line) or t' = 0.5 (dashed line), (b) BCC lattice 
with t' = (solid line) or t' = 1.0 (dashed line), and (c) SC 
lattice. 



In introducing the interaction, we start with a discus- 
sion on how strong the Coulomb interaction may be for 
the FLEX approximation to be reliable. In Fig. 3, we 
plot the static magnetic susceptibility, x = dTr[a z G]/dh, 
and the RPA spin susceptibility, xrpa(0,0), as a func- 
tion of temperature T for t' — 0.5, n = 0.2 for various 
values of U = 0.4 - 2.0. We deliberately took t' = 0.5 
and n = 0.2 because the ferromagnetic spin fluctuation 
turns out to be the strongest for this parameter set as 
we shall see below, so that FLEX is expected to be least 
reliable. Indeed, the difference between \ an d Xrpa is 
seen to be large for [7 = 2. This means that we cannot 
approximate the irreducible effective interaction or the 
four-point vertex T = 5 2 $>/5G5G as U in this case. 

As a further check, Fig. 4 looks at the double occu- 
pancy, (n^ni) normalized by (r7,|}(n^), evaluated from 
the relation, TrG a Ilc, = U(n^ni). We can see that for 
U = 2 and n = 0.2, (nfnj,) becomes negative, which is 
unphysical. Hence, for t' = 0.5 and n = 0.2, a most 
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Fig. 3. The FLEX static magnetic susceptibility (dashed lines) 
and the RPA susceptibility (solid lines) for the Hubbard model 
on an FCC lattice with t' = 0.5 and n = 0.2 as a function of 
temperature for U = 2(a), U = 1(b) and U = 0.4(c). 



stringent case, the FLEX approximation breaks down at 
least for U > 2. Thus we set U < 2 hereafter. 

Fig. 4. A reliability check for the FLEX method from the double 
occupancy (rif rij) normalized by (n-r){ri{) as a function of U for 
t' = 0.5 FCC lattice for various values of the band filling, n. 



Now we are in position to present the result for the 
density of states and susceptibility. Figure 5 shows the 
DOS for t' = 0.5, U = 1,2 for the density varied from 
n = 0.2 (a), n = 0.6(b) to n = 1.0(c) with T = 0.1. The 
Fermi level is seen to be right at the peak of the DOS 
for n = 0.2, so that the ferromagnetic state should be 
favored according to Stoner's mean field picture. 

Fig. 5. The density of states (integrated one-particle spectral 
function) for t' = 0.5 FCC lattice with n =0.2 (a), 0.6(b), 1.0(c). 
The solid (dashed) lines are for U = 1(2). The insets are blow- 
ups of the DOS around the Fermi level. 



If we then look into the RPA FLEX spin susceptibil- 
ity XRpA(k, 0) in Fig. 6, the ferromagnetic spin fluctu- 
ation is indeed strong (i.e., x(k) enhanced at T point) 
for n — 0.2 as expected from the DOS. (Note that the 
susceptibility takes the same value at T, k = (0, 0, 0) and 
K, (it, 7t, it) for the convention of taking the FCC disper- 
sion (2).) 

Fig. 6. The RPA spin susceptibility as a function of the 
wavenumber over the Brillouin zone (inset) for t' = 0.5 FCC 
lattice with n =0.2 (a), 0.6(b) or 1.0(c) for U = (solid lines), 
U = 1 (dashed lines) or U = 2 (dotted lines). 



For an intermediate density n — 0.6, on the other 
hand, the ferromagnetic spin fluctuation .-disappears. 
This is in contrast with the results by Ulmke,^ who pre- 
dicted, for the strong-coupling regime, that the ferromag- 
netic state is most favored for intermediate densities, and 
that there is a finite lower critical density n c = 0.1 ~ 0.2 
below which the ferromagnetic state is not realized. For 
n = 1.0 (half filling), we cannot find any instabilities at 
(tt, 0, 0) either, which has been predicted by Ulmke in 
the strong-coupling regime S3 

In Fig. 7, we show xrpa(0, 0) for U = 2 as a function 
of n. We can see that the ferromagnetic fluctuation is 
most dominant for n = 0.2 ~ 0.3. If we calculate, for 
U = 5, the temperature at which xrpa(0, 0) diverges as 
a function of the density n, we find that the temperature 
takes its maximum value ~ O.lt for n ~ 0.2. Since the 
value of U — 5 is too large for the FLEX approximation 



to be valid, at least for n ~ 0.2 as mentioned above, 
we cannot naively identify these as the transition tem- 
peratures to a ferromagnetic state. However, the den- 
sity dependence of these temperatures suggests that the 
ferromagnetic spin fluctuation is most dominant around 
n ~ 0.2. Therefore, we believe that the tendency of 
the ferromagnetic spin fluctuation being most dominant 
around n ~ 0.2 in the weak coupling regime should be 
real. 

Fig. 7. XRPa(0) plotted as a function of n for (7 = 2 and T = 0.1. 



Having identified in FLEX that the ferromagnetic spin 
fluctuation is the strongest for FCC with n ~ 0.2, we 
move on to check whether the strong ferromagnetic spin 
fluctuation remains in the TPSC approximation, which 
includes vertex corrections. The qualitative features of 
the spin susceptibility x obtained with TPSC has turned 
out to be similar to Xkpa obtained with FLEX. Here, we 
focus on how strong the ferromagnetic spin fluctuation 
diverges for the case of n = 0.2 and t' = 0.5. We have 
checked the internal self accuracy-ie., the difference be- 
tween Tr£«G*W and TrSWG^P to find around 10% 
accuracy at worst. To check the convergence of eq. (19) 
or eq. (20), 8192 Matsubara frequencies are required, so 
we have calculated up to N = 16 3 . 

The TPSC result for 1/x as a function of T is plotted in 
Fig. 8 for U = 1 ~ 3. Here the chemical potential /i is ap- 
proximated to fiQ. The result exhibits a Curie- Weiss like 
behavior ( X ~ 1/(T + 9)) with 9 > for U < 3, so that 
the ferromagnetic transition does not take place as far 
as this result is concerned. If we approximate the chem- 
ical potential /i as fi\ — E(kp,0), as mentioned at the 
bottom of §IIC, we can see that the ferromagnetic spin 
fluctuation becomes more divergent (6 < for {7 = 3) 
as seen in the inset. For t' = 0.5 and n — 0.2, a small 
change in the chemical potential causes large changes in 
the solution for U sp or C/ C h (eqns(19,20)). However, the 
tendency detected here that the chemical potential de- 
creases as the approximation for the chemical potential is 
improved should be correct, and it can be shown, at least 
in the present case of FCC, that U sp becomes larger for 
smaller fi. Therefore, although /i = \i\ — £(kp, 0) suffers 
from some numerical inaccuracies due to an error in the 
Pade fit to £(kp, 0), the dominance of the ferromagnetic 
spin fluctuation should be reliable. 

Fig. 8. The inverse of spin susceptibility for FCC in the TPSC 
approximation as a function of T for U = 1(d), U = 2(+) or 
U = 3(x), where [iq is adopted for the chemical potential. The 
inset is a similar plot, where /ii — S(kj?,0) is adopted for the 
chemical potential 



When we turn off the next nearest neighbor hopping 
(f = 0), the DOS's obtained by FLEX for n = 0.2(a), 
n = 0.6(b), and n — 1.0(c) in Fig. 9 show that the peak 
around the bottom becomes much weaker than that for 
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t' ^ 0. Accordingly, XRPA(k, 0) obtained with FLEX in 
Fig. 10 shows that the ferromagnetic spin fluctuation 
at low density (n ~ 0.2) is much weaker than that for 
t' = 0.5. We can in fact quantify this from the fact that 
only when U is taken to be as large as 10 (where FLEX 
is not reliable) does 1/xrpa(0, 0) become zero for t' = 0. 

Fig. 9. The density of states for t' = FCC lattice with n =0.2 
(a), 0.6(b), 1.0(c). The solid (dashed) lines are for U = 1(2). 
The insets are blow-ups of the DOS around the Fermi level. 



Fig. 10. The RPA spin susceptibility as a function of the 
wavenumber for t' = FCC lattice with n =0.2 (a), 0.6(b) or 
1.0(c) for U = (solid lines), U = 1 (dashed lines) or U = 2 
(dotted lines). 



3.2 BCC lattice 

We next turn to the BCC lattice, a typical bipartite 
lattice. We again study the two cases in the presence and 
absence of the next nearest neighbor hopping, t'; (i) a fi- 
nite t' = 1.0 with n = 0.3, 0.7, 1.0, and (ii) nearest neigh- 
bor hopping only (f = 0.0) with n = 0.2, 0.6, 0.8, 1.0. We 
set T = 0.1 as before. First, the non-interacting DOS, 
displayed in Fig. 2(b), has a peak around the band cen- 
ter, which resides exactly at the center for t' = 0.0, since 
a bipartite lattice has an electron-hole symmetry when 
electrons hop to nearest neighbors only. As for the accu- 
racy check, ("t rl l)/( n T}( n |} remains positive for U < 2 
for cases (i) and (ii). 

We start with the result for a finite t 1 = 1.0. If we look 
at the DOS for n = 0.3(a), n = 0.7(b), and n = 1.0(c) in 
Fig. 11, for n ~ 0.7 the Fermi level is located around the 
peak of the DOS, although a fictitious dip appears due 
to a numerical inaccuracy in the analytic continuation. 

Fig. 11. The density of states for t' = 1.0 BCC lattice with 
n =0.3 (a), 0.7(b), 1.0(c). The solid (dashed) lines are for 
U = 1(2). The insets are blow-ups of the DOS around the Fermi 
level. 



In Fig. 12, we plot xrpa obtained with FLEX for 
n = 0.3(a), n — 0.7(b), and n = 1.0(c). Although a finite 
size effect (slight wiggles) is visible even when we take 
N = 32 3 , we believe that the overall feature of the spin 
structure is reliable. As expected, the ferromagnetic spin 
fluctuation is dominant when Ep resides at the peak for 
n ~ 0.7, but the enhancement in \ is much weaker than 
that for the FCC lattice with t' = 0.5 and n = 0.2.0 
In fact 1/xrpa(0, 0) does not become zero for any value 
of U for T = 0.1. Although the divergence of the peak 
in the DOS is as sharp as that of FCC with t' = 0.5, 
the strength of the ferromagnetic spin fluctuation differs 
between these cases. This implies that the peak of the 
DOS located around the bottom is more advantageous 



for the ferromagnetic spin fluctuation than when located 
around the center. 

Fig. 12. The RPA spin susceptibility as a function of the 
wavenumber for t' = 1 BCC lattice with n =0.3 (a), 0.7(b) 
or 1.0(c) for (7 = (solid lines), U = 1 (dashed lines) or U = 2 
(dotted lines). 



When we turn off t' in the BCC lattice, the Fermi 
level coincides with the peak of the DOS for n = 1.0 as 
shown in Fig. 13. Here again, the analytical continuation 
with the Pade approximation is unstable for n > 0.6, 
causing a dip to appear in the DOS. In Fig. 14, we 
plot the xrpa obtained with FLEX for n = 0.2 ~ 1.0. 
This time the spin fluctuation becomes dominant at the 
artMerromagnetic point (X point with k = (tt, 7t, it)) for 
n = 0.8 and n — 1.0. For n — 1.0, 1/xrpa(0, 0) becomes 
zero for U ^2 and T = 0.1. As we shall see below, the 
antiferromagnetic fluctuation is much stronger than that 
for the simple cubic lattice. 

Fig. 13. The density of states for t' = BCC lattice with n =0.2 
(a), 0.6(b), 0.8(c) or 1.0(d). The solid (dashed) lines are for 
U = l{2). 



Fig. 14. The RPA spin susceptibility as a function of the 
wavenumber for t' = BCC lattice with n =0.2 (a), 0.6(b), 
0.8(c) or 1.0(d) for U = (solid lines), (7 = 1 (dashed lines) or 
U = 2 (dotted lines). 



3.3 Simple cubic lattice 

We finally turn to the SC lattice. We plot the DOS 
in Fig. 15 for n = 0.2, n = 0.6, n = 1.0, while the DOS 
for the non-interacting case is shown in Fig. 2(c). They 
are electron-hole symmetric since SC lattice is bipartite. 
We set T = 0.1. Here again (^t rl i)/( rl T)( n l) remains 
positive for U < 2. If we look at xrpa obtained with 
FLEX in Fig. 16, the ferromagnetic fluctuation is seen 
to be suppressed. The only enhancement appears for 
n = 1.0 at the antiferromagnetic point (K point), but 
the fluctuation is weaker than that for the BCC lattice 
above. 

Fig. 15. The density of states for t' = SC lattice with n =0.2 
(a), 0.6(b), or 1.0(c). The solid (dashed) lines are for U = 1(2). 



§4. Discussions and Summary 

We have studied the Hubbard model on 3D FCC, 
BCC, and SC lattices with the FLEX and the TPSC 
approximations in a weak-coupling regime. We have cal- 
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Fig. 16. The RPA spin susceptibility as a function of the 
wavenumber for t' = SC lattice with n =0.2 (a), 0.6(b) or 
1.0(c) for U = (solid lines), (7=1 (dashed lines) or U = 2 
(dotted lines). 



culated the spin susceptibility as well as the density of 
states (integrated single-particle spectral function). For 
cases in which the DOS is divergent near the bottom 
of the band and the Fermi level is located around the 
divergence, the ferromagnetic instability is found to be 
dominant. Such a situation is realized for a low band fill- 
ing in the FCC lattice with a finite next-nearest neighbor 
hopping (t 1 ) or close to the half-filling in the BCC lattice 
with a finite t' . For bipartite lattices (SC and BCC lat- 
tices with nearest neighbor hopping only), by contrast, 
the antiferromagnetic fluctuation is dominant around the 
half-filling. 

We have found, among the cases studied, that the 
strongest ferromagnetic fluctuation occurs for a low band 
filling FCC lattice with t' , whose DOS is characterized 
by the sharpest divergence that resides around the band 
bottom. From a view of "frustration" we may say the fol- 
lowing. In the absence of t' , FCC is frustrated in that its 
constituent plaquette is a tetrahedron where antiferro- 
magnetic correlations have to interfere, while BCC is bi- 
partite compatible with antiferromagnetism. Thus they 
have different situations, before t is turned on to make 
the DOS for BCC divergent. While both of the diver- 
gent DOS and the frustration can favor ferromagnetism 
intuitively, we cannot at present stage distinguish effects 
of the two. We can still observe that the present result 
for a stronger ferromagnetic fluctuation in FCC than in 
BCC implies that the peak of the DOS lying around the 
bottom of the band may be a key factor. As stressed 
in Introduction, such a comparison could not have been 
performed from Kanamori's theory, which ipcuses on di- 
lute cases. For ground state, Hanisch et a/.EP suggested 
that the fully-polarized ground state is favored when the 
DOS diverges near the bottom of the band. On the other 
hand, we have discussed here the tendency toward ferro- 
magnetic instabilities at the finite temperatures, which 
includes the possibilities of not only fully-polarized state 
but also partially polarized state. Hence the relation be- 
tween the results of Hanisch et al. and our results is not 
trivial and is an interesting future problem. 

The discussion on the position of the divergence in 
DOS and frustration reminds us of a situation in ID, 
in a quite different context. In the flat-band ferromag- 
netism mentioned in Introduction, we can consider quasi- 
one-dimensional systems in which flat bands are real- 
ized. There, flat bands can arise in non-frustrated models 
(such as a chain of squares) but the flat band resides 
in the middle of the band, so that the ferromagnetism 
is shown to be destroyed when U becomes too large. In 
frustrated models (such as a chain of triangles), by con- 
trast, the flat band is situated at the bottom, and the 
ferromagnetism is guaranteed for < U < oo from the 
Mielke-Tasaki theorem. Although this observation con- 
cerns quite a distinct class of systems, further discussions 



may prove fruitful. 

The FLEX result for the FCC lattice indeed shows that 
the ferromagnetic instability is most dominant in the 
low-density region. This reminds us that Ni has an FCC 
structure with a low hole concentration on the one hand. 
On the other, the result is contrasted with a DMFT re- 
sult for the strong-coupling regime by UlmkeJiS' which 
suggests that for FCC lattice with t' ^ 0, the ferromag- 
netic region only appears above n c — 0.1 ~ 0.2 and 
is peaked at an intermediate density (n ~ 0.6). It is 
an interesting future problem to explore how these two 
regimes are interpolated. 

The antiferromagnetic fluctuation, by contrast, ap- 
pears for bipartite (i.e., electron- hole symmetric) lattices 
(BCC or SC), and the fluctuation is the strongest for the 
half-filled BCC with a divergent DOS at the Fermi level 
(for t' = 0) of all the situations surveyed here. We note 
that it is stronger than that in the half-filled SC lattice 
whose DOS has only a finite peak. This cannot be un- 
derstood from x p h(7T, 7r, 7r) for non-interacting electrons, 
since there is not much difference in this quantity be- 
tween the two cases. 
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